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ABSTRACT 

We calculate the Isgur-Wise function from the solutions of the Bethe-Salpeter equa- 
tions. The shape of the Isgur-Wise function thus calculated is a prediction of the 
Bethe-Salpeter equations and does not depend on undetermined parameters. We de- 
velop an analytical approximation to our Isgur-Wise function in the form ^{uj) = 
r][l - ^{u - 1) + a{uj - 1)3/2] ^j^g^g p2 ^ 1279, a = .91, r] = .9942 and u is the 
recoil velocity. The Isgur-Wise function is then used to obtain Vcb from the recent exper- 
imental data of B D*iu decay. Our best estimate of Vcb is (34.7 ± 2.5) x 10^^, which 
is comparable to some of the latest estimates in the literature. 



1. INTRODUCTION 



Heavy quark effective theory (HQET) |[T|, H, ^ |] has opened a new window for the 
determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [§, ^, ^. In 
particular, determination of Vcb using theoretical predictions of HQET and experimental 
measurement of the differential decay rate of the exclusive semileptonic decay B D*iu 
has received a great deal of attention. With the recent high-precision data available from 
several experimental groups 0, ||, Q , it is now possible to determine the value of Vcb with 
reasonable accuracy. 

In order to extract Vcb from the experimental data, one needs to calculate hadronic 
form factors which include nonperturbative effects. HQET has vastly simplified these 
calculations and only one universal function, called the Isgur-Wise (IW) function 
^, |10|, has been shown to play a central role in many calculations involving decay of 
heavy mesons. The calculation of the IW function however is model dependent and 
several different parameterizations for it have been used in the literature 0, P). After a 
particular parameterization is chosen, one typically fits the experimental data to extract 
the unknown parameter (s) in the IW function as well as Vcb from the experimental data. 

In this paper we use a covariant reduction of the Bethe-Salpeter equation (BSE) 
to calculate the IW function |[12|. The BSE was solved numerically and the parame- 
ters appearing in it (namely the quark masses, string tension and the running coupling 
strength for one-gluon exchange) were determined by fitting the calculated spectrum to 
the observed masses of more than 40 mesons []T3|, The resulting meson mass spec- 
trum agrees very well with the experimental data. Once the parameters are thus fixed, 
the meson wave functions from the BSE can be used to predict physical observables. In 
particular, the IW function may then be evaluated from the wavefunctions of the BSE 
and would represent a prediction independent of any undetermined parameters. Knowing 
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the shape of the IW function we can then extract Vcb from the experimental data 

This paper is organized as follows: in Section 2 we recall the main ingredients for the 
calculation coming from HQET. Section 3 describes the salient features of the BSE and 
describes the meson wavefunctions of interest. In Section 4 we calculate the IW function 
and extract Vcb from the experimental data. We conclude the paper in Section 5 with a 
brief outline of future work. 

2. HEAVY MESON DECAY FORM FACTORS 



The hadronic matrix elements for the decay B ^ D* £u take a simple form when de- 
scribed in the context of HQET. Such decays are mediated by heavy quark currents 
Vfj. = c'-f^b and = cr^^'^^h and the corresponding matrix elements are in general de- 
scribed in terms of four form factors denoted by i^j, i = 1, 2, 3, 4 : 

{D*{v',t)\V^'\B{v)) = zy^^^^7^ei(c^)e^""^e>>/3 

{D\v\e)\A^\B{v)) = ^rriBmnM^) {uj + 1) e*^-{U^) + U^)v'^y.v, (2.1) 

where uj = v.v', v and v' being the velocities of B and D* meson respectively. 

In the limit where masses of the heavy quarks tend to infinity, the form factors 
satisfy the conditions 

ei = 6 = e4 = eM'), 6 = 0, (2.2) 

where ^{v.v') is a single universal function, called the IW function [0, 0. In the limit of 
infinitely heavy quark masses the IW function is normalized to unity at zero recoil, i.e. 

e(i) = i. 

The IW function can be related to the overlap integral of normalized meson wave 
functions in the infinite momentum frame. If ipis and ipin* denote the wavefunctions of 
the light degrees of freedom in B and D* mesons respectively, then the IW function can 
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be written as 



= (^)'^V i^!D^iv')i^iBiv)d'x. (2.3) 



In the heavy quark hmit close to u; = 1, the IW function has the form 

eM = l-p2(c.-l) + 0[(c.-l)2], (2.4) 

where is the slope of the Isgur-Wise function at = 1. 

The differential decay rate for the process discussed above is given by p[ 
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^{rriB — mD* Y m\), rj\ \/uS^ — \ {uo + 1)^ 



da; 487r3 

4c<j 1 — lujr + 



X 



1 + 



K^rei^), (2.5) 



+ 1 (l-r)2 

where r = and ?7a is a constant which is present due to a finite renormalization of 
the axial vector current. 

It is clear from the above expression that the knowledge of the differential decay rate 
and of the Isgur-Wise function would allow us to calculate Vch- However, as is evident 
from Eqn.(2.3), we need to know the meson wave functions to calculate the Isgur-Wise 
function. In the next section we show how to do this using a covariant reduction of the 
BSE. 

3. MESON WAVEFUNCTIONS FROM THE BSE 



In a previous set of works, reductions of the BSE has been solved for qq systems 



13, 14|. In this section we shall briefly describe our treatment and the wavefunctions 



resulting from our calculation. The actual equations solved to obtain these wavefunctions 
are given in the appendix and the reader is referred to the previous works cited above for 
a full treatment. 

The BSE is a covariant four-dimensional wave equation for relativistic bound states 
and is very challenging to solve exactly for realistic kernels. One typically uses several 
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approximations to reduce it to a solvable form. We have used a quasipotential equation 
framework to reduce the BSE to the three dimensional integral equation introduced in 
[|T3| and have retained only the ladder diagram component of the full BSE kernel. Under 
this approximation the interaction kernel contains a one-gluon exchange term Voge, to 
which we add a phenomenological, long-range confining potential, Vcon- The interaction 
kernel that we use is thus derived from the potential 



+ v^on = r^j^^. + -]!i^odjr^^(^^^iyT7- (3-1) 



Here, is the strong coupling, which is weighted by the meson color factor of |, and the 
string tension a is the strength of the confining part of the interaction. 

We also take into account the running of the strong coupling constant. Specifically, 
we use the form 

"^^"^^ 4vr + /3a.(/i2)ln(gV/x2)' ^^'^^ 

where 

Q' = 7'Ml_ + f3', (3.3) 

and where 7 and (3 are parameters determined by a fit to the meson spectrum. In the 
above equation fi is taken to be the mass of Z boson and Q;s(/i^) is correspondingly chosen 
to be equal to 0.12 based on experimental measurements. 

In our formulation of BSE there are seven parameters : four masses, mu=Tnd, rric, 
rris, mh] the string tension a, and the parameters 7 and (3 used to govern the running of 
the coupling constant. These parameters were determined by fitting the meson masses 
calculated from the BSE to the observed spectrum. We utilize this BSE model for the 
mesons to evaluate the meson wavef unctions. The wavefunctions for B and D* mesons 
are shown in Figs. (1) and (2) respectively. For the purpose of further calculations we 
have also obtained analytic representations of these wavefunctions : 

^(x,,,.) = ((^^)^e--(^W)^, (3.4) 



where, A = .59 GeV, n = 1.16 for the B meson and A = .42 GeV, n = 1.52 for the 
D* meson. We have plotted the numerical wavefunctions and their analytic representa- 
tions together to display the accuracy of the latter. For convenience, we shall use these 
analytical expressions for subsequent calculations in this paper. 

4. ISGUR-WISE FUNCTION AND 



In this section we shall describe our results using the ingredients that have been 
presented above. Using Eqn.(2.3) and taking into account the relativistic boost (assumed 
along the z direction), the IW function can be written as 

e(^) = i^y^' f ^tn,{x,y,u;z)^i^{x,y,z)e'''^^^'d\ (4.1) 
iO + 1 J 

where E is the mass of the light component of the D* meson in the rest frame of the B 
meson. 

We use the wavefunctions derived from the BSE to evaluate the IW function and we 
emphasize again that our IW function involves no additional parameter fitting. The plot 
of our IW function is shown in Fig. (3). It is interesting to note that since the wave 
functions used are coming from the BSE which is solved independent of the heavy quark 
approximation, the IW function does not go to unity at zero recoil. The deviation of the 
IW function from unity at zero recoil can be attributed to the finite mass corrections which 
are incorporated in our BSE model. We have also obtained an analytical representation 
of the IW function that is plotted in Fig. (3). This is given by 

e(^) = V[l - ^(^ - 1) + a{u; - 1)3/2], ^4 2) 

where = 1.279, a = .91 and rj = .9942. 

We can now use the IW function as calculated above to extract Vcb from the experi- 
mental data ioi B ^ DHv decay. We have used the data from ARGUS 93 0, CLEO 93 
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Table 1: Values of Vcb in the units of 10~^ from different sets of experimental data. The 
first row indicates the values of Vcb obtained by using the IW function that was calculated 
numerically. The second row indicates the corresponding values of Vcb obtained by using 
an analytical approximation to the IW function as given in Eqn. (4.2). 





ARGUS 93 


CLEO 93 


CLEO 94 


From Numerical IW Function 


(35.8 ±8.3) 


(34.8 ±6.1) 


(34.7 ±2.5) 


From Analytical IW Function 


(35.9 ±8.3) 


(34.8 ±6.1) 


(34.7 ±2.5) 



P and CLEO 94 |^ and the corresponding fits to these data using our IW function are 
shown in the Figs. (4), (5) and (6) respectively. 

In Table 1 we present the different values of Vcb (in the units of 10^'^) obtained from 
these experimental data (for each set of experimental data we have shown two values of 
Vcb, one obtained from the numerical calculation of the IW function and the other from 
its analytic representation as given in Eqn. (4.2).) The uncertainties in the values of 
Vcb tabulated above are determined by fitting to the two extreme values appearing on the 
experimental error bars. The different values thus obtained are consistent with each other 
and are comparable to some of the latest estimates of Vcb B H; il • 

5. CONCLUSION 

We have used the solutions of the BSE to predict the shape of the IW function. Such a 
prediction is independent of any undetermined parameter. The solution of BSE does not 
depend on the heavy quark approximation and the corresponding IW function calculated 
from its solutions incorporates some of the finite mass corrections. This is evident from 
the fact that in our calculation the IW function deviates from unity at zero recoil. The 
analytical representation of the IW function shows an interesting power law behaviour 
which is different from the usual representations available in the literature. Using the 
IW function we have calculated Vcb from the latest set of experimental data. Our best 
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estimate of Vcb = (34.7 ±2.5) x 10^'^ is comparable to the latest estimates of Vcb available 
in the literature ||, |[. 

From the knowledge of the solutions of Vcb it should be possible to calculate the 
four form factors in Eqn. (2.1) directly without using the heavy quark limit. Such a 
calculation would provide an interesting test of the accuracy and consistency of the heavy 
quark approximation and is currently being pursued. 
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APPENDIX 



A. BETHE SALPETER EQUATIONS 

In this appendix we present the BSE that have been solved. This is done for the sake 



of completeness and the reader should refer to [|T^, [T^ for further details. 
The BSE for the 0~ channel is given by 

[{E, + E,f - E']i,\q) = -[ir^%q') + /o^V°(g')], (A.l) 

vrg 

and those for the 1~ channel are given by 

[{E, + E,f - E']i,'\q) = ^[{ir + /r)^^'(g') + (/r + l2'V{q')], (A.2) 



and 



[{E, + E,y - E']r\p) = -[{ir + ir)i^'\q') + {ir + nm^'^Q)]- (a.3) 



The symbols appearing in the above equations are defined as : 

/r^^lim^ dq(l[Qo{z')Tr + Ql{z)Tri 
ir = <^lini— / dq q [QUz )Tr], 
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ir = ^lini— / cig'g'[gi(^')Tr + Ql{z)Tr], 

IJ.-*0OjJ,'^ Jo 

ir ^^-^ r dqq[Qo{z)Tr + Qliz)Tr], 
6 JO 

ir r dqqmz)Tr + Qi{z)Tr], 

6 Jo 

ir=^-^ rdqq[Ql{z)Tr], 
6 Jo 

ir-^-^ rdqq[Ql{z)Tr], 
6 Jo 

ir r dqqmz)Tr + Qi{z)Tr], 

6 Jo 



where 



Tr = j^^^^^^^{E,E,E[E', - q'q' + g^i^l^; - q^E.E,), 

T"s^ ^ ^ , , nq'(Ei - E2)(E[ - E'), 

^ 4:{EiE2E[E2) ^ ^ ^ 

Tr = j^^^^^^^{E,E2E[E', - q\" - q^E'X + q'e.E^), (A.4) 



E,^{mi + q')-^, E2^{mi + q')-^, 

0/2 9/2 9 

? + ? / g"' + ? + //^ 

= — 1—. z 



2qq 2qq 
In the above equations, Qj are the Legendre polynomials of the second kind and 

Q]{z) = zQj{z) - 6jo, 

Q'(z)^-^{jzQ,{z) + Q,.,{z)), 



J + 

Q'jiz)^i-^)HzQjiz) + Qj^,{z)). 



Finally, we also have 

^{-J + 1) 



= 'V + ^X?) + J^=^' ^■(?)). (A.5) 

(2(2j + l))2 

i/j^'^^iq) and '^^^{q) are obtained by numerically solving Eqns. (A. 2) and (A. 3). 
iljj^j_i{q) and iljj^j_^_i{q) are then obtained numerically from Eqn. (A.5). These are the 
wavefunctions for the D* meson. The / = 2 component of the D* meson wavefunction 
is negligible in comparison to the I — component and is neglected. Similarly, ip'^{q) 
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(representing the wavefunction of the B meson) is obtained by numerically solving Eqn. 
(A.l). Analytical approximations of ip^{q) and ipl{q) are given in Eqn. (3.4). They are 
used in Eqn. (4.1) to calculate the IW function. Fig. 1 and Fig. 2 show the numerical 
wavefunctions and their analytical representations for B and D* mesons respectively. 
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FIGURE CAPTIONS 



Figure 1: Radial wavefunction for tlie B meson. Tlie solid line represents the numerical 
solution of the BSE. This is obtained from the Fourier transform of ip^iq) appearing in 
Eqn. (A.l). The entire Fourier transform of ip^iq) in position space is normalized to 
unity. An analytical approximation to ip^iq) is given in Eqn. (3.4) (with A = 0.59 GeV 
and n = 1.16) and is represented by the dashed line. 

Figure 2: Radial wavefunction for the D* meson. The solid line represents the numerical 
solution of the / = component. This is obtained from the Fourier transform of ipl{q) 
appearing in Eqn. (A. 5). The entire Fourier transform of ipoi^) position space is 
normalized to unity. An analytical approximation to ipoi^) given in Eqn. (3.4) (with 
A = 0.42 GeV and n = 1.52) and is represented by the dashed line. 

Figure 3: The Isgur-Wise function. The solid line represents the IW function calculated 
from Eqn. (4.1). An analytical approximation to the IW function is given in Eqn. (4.2) 
and is represented by the dashed line. The deviation of the IW function from unity at 
zero recoil is due to the finite heavy quark mass effects incorporated in our BSE model. 
Figure 4: Plot of Kfc ^{uj) vs. u from ARGUS 93 data of B ^ D*iu decay. 0. The 
data points for Vcb ^{uj) are extracted from the experimental result by using Eqn. (2.5). 
The solid line represents our best fit to the data. 

Figure 5: Plot of ^{oo) vs. u from CLEO 93 data of S ^ D*iiy decay. 0. The data 
points for Vcb ^{^) are extracted from the experimental result by using Eqn. (2.5). The 
solid line represents our best fit to the data. 

Figure 6: Plot of ^{oo) vs. uj from CLEO 94 data of B ^ D*iiy decay. 0. The data 
points for Vcb ^{^) a^re extracted from the experimental result by using Eqn. (2.5). The 
solid line represents our best fit to the data. 
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